function radius = calculateRadius(Vl)
%CALCULATERADIUS A function used in substationlocation.m

sizetmp=size(Vl);
length=0;
for i=1:sizetmp(1)
	for j=i+1:sizetmp(1)
		lengthtmp = (Vl(i,1)-Vl(j,1))^2+(Vl(i,2)-Vl(j,2))^2;
		if lengthtmp>length
			length=lengthtmp;
		end
	end
end
radius = sqrt(length)/2
